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Abstract 

Charm production as related to the determination of the strange sea density in deep in- 
elastic charged current processes is studied predominantly in the framework of the MS 
fixed flavor factorization scheme. Perturbative stability within this formalism is demon- 
strated. The compatibility of recent next-to-leading order strange quark distributions 
with the available dimuon and F£ N data is investigated. It is shown that final conclu- 
sions concerning these distributions afford further analyses of presently available and/or 
forthcoming neutrino data. 



Heavy quark production at high energy neutral current reactions was recently shown 
to be optimally described within the framework of a fixed flavor (factorization) scheme 
(FFS) where, besides the gluon g, only the u, d and s quarks are considered as partons 
and any heavy quark (c, b, ...) contribution is calculated in fixed order a s perturbation 
theory. Within this framework the charged current production of a heavy quark pair 
such as tb in up —> fi~ tb X follows the same pattern || utilizing the relevant formulae 
for the underlying Wg fusion' subprocess W + g — > tb. Since both m tt b 3> Aqcd, we do 
not encounter any mass singularities here and a treatment within the framework of the 
FFS is straightforward and unproblematic. This favorable situation changes, however, 
when considering the corresponding charm production process (e.g. W + g — ► cs) since 
the associated strange quark is taken as massless in the FFS and considered as a parton. 
In contrast to the former cases we encounter here a mixed situation which affords a 
careful treatment within the framework of the FFS. The leading order (LO) contribution 
for charm production in uN — > /x~ c X, N = (p + n)/2, comes from the basic O{a s ) 
subprocess W + s' — »■ c where 
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8 V N = \ V cs\ S + \V cd \ — — (1) 

with |V^ S | = 0.9743 and \V c d\ = 0.221. The W + g —> cs fusion process yields the essential 
part of the next-to-leading order (NLO) correction where the other part is due to the sub- 
process W + s' — > gc. Both subprocesses posses a mass singularity associated with m s = 
which is absorbed via dimensional regularization into the renormalized, Q 2 -dependent, 
parton distribution s'. The remaining finite pieces of W + g — > cs and W + s' — > gc then 
yield the genuine NLO correction to W + s' — > c which will henceforth be considered in 
the now commonly adopted MS factorization scheme. Denoting the contributions of the 
above subprocesses to the structure functions Fi(x, Q 2 ) by F?(x, Q 2 ) and defining further- 
more T{ = Ff, Tl = F 3 c /2, Tl = F^jli where £ = x(l + m 2 JQ 2 ), one obtains in NLO 
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Here A = Q 2 /{Q 2 + m 2 ) and the Hf 9 are given, up to minor modifications specified in 
the Appendix, in ref. ||. The specific choice for the factorization scale /i will be studied 
below. The inclusive cross section in terms of the Fi(x, Q 2 ) is given by 
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To study the size of the NLO corrections we shall utilize the LO and NLO parton 
distributions of Q which are already conceived in the FFS being furthermore MS distri- 
butions in the NLO. In addition we shall also employ the LO and NLO(MS) CTEQ3 |J 
and the NLO (MS) MRS (A) || parton densities which refer to the 'variable flavor' scheme 
where intrinsic charm densities are purely radiatively generated using the ordinary mass- 
less evolution equations, starting at Q — m c . For definiteness we show in fig. 1 the 
quantity 

* 8 & Q )e// = 2 G\M N E V ^ Ix-Jy- (4) 
which has been also studied experimentally |?| and where the superscript cs refers just to 
the CKM non-suppressed (V cs ) component of s' in eqs. (1-3). Note that in LO the cross 
section in (4) reduces to 
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As can be seen in fig. 1 the NLO corrections to the LO results are reasonably small and, 
in particular, do not afford a drastic change of s(x, Q 2 ) when passing from the LO to the 
NLO analysis. 

This contrasts with the conclusions of the CCFR group whose NLO s(x, Q 2 ) is 
almost twice as large as compared to their previous || LO s(x,Q 2 ). The analysis of the 
CCFR group is based on the NLO formalism of || which is not strictly equivalent to 
our NLO (MS) FFS formalism but still is expected to yield quite similar results. The 
enhancement of the NLO s(x, Q 2 ) in can therefore not be attributed to the different 
formalism itself [T(J but rather to its inconsistent application. In the formalism of ref. 



one considers the W + g — > cs contribution with m s ^ 0, i. e. employs a finite mass 



regularization and subtracts from it that (collinear) part which is already contained in 
the renormalized, Q 2 -dependent s(x,Q 2 ). The CCFR group applied their acceptance 
corrections to the full contribution from W + g — > cs, which corresponds effectively to a 
multiplication with an acceptance correction factor A = 0.6 ±0.1, while inconsistently 
keeping the subtraction term in its full original (acceptance uncorrected) strength fTT| . 
This latter subtraction term is given, relative to £s(£, n 2 ), by 

«m to 4 r* S iz,s)pM(i) (5) 

27T mi h z qy \z 



ft \ i ' ] 

SUB 



using |l Ij m s = 0.2 GeV. In fig. 2 we compare the result obtained in this manner with 
the one where also the subtraction term || in (5) was consistently multiplied by the same 
acceptance factor A. The result corresponding to the acceptance uncorrected subtraction 
term (dashed curve) clearly demonstrates that SUB alone [eq. (5)] represents too strong 
a suppression of the mass singularity component in LO + A * NLO and that the correct 
result (solid curve) in fig. 2 is almost a factor of 2 larger in the small- x region. This 
implies that instead of snlo — 2 slo for x ~ 0.1, as inferred by CCFR |7|] and used 
for our analysis in fig. 2, one rather needs a smaller snlo> i- e - closer in size to slo> i n 
order to reduce the solid curve in fig. 2 and to bring it closer to experiment. Here we 
have chosen @ a factorization scale fi = 2 p™ ax = A(W 2 , m 2 c , M%)/W, where p™ ax is the 
maximum available transverse momentum of the final state charm quark; the results in 
fig. 2 remain practically unaltered with the alternative choice fi 2 = Q 2 + rr? c . 

A further feature emanating from the fits in was m^ LO ~ 1.7 GeV as compared 
to the previous m^° ~ 1.3 GeV which further suppressed the NLO cross section and 
demanded the unusual, even more enhanced NLO s(x, Q 2 )- A consistent treatment of the 
acceptance correction would most probably also lower the fitted m^ LO down to a more 
reasonable m^ LO ~ 1.5 GeV and bring the NLO s(x, Q 2 ) close to the LO s(x, Q 2 )- 

Our conclusions concerning the strange quark distributions of j|, [I| are that they 
agree in LO with || and are not refuted in NLO by the analysis in [0]. Furthermore 
due to the perturbative stability demonstrated in fig. 1 we expect the NLO strange quark 
distributions of [4-6] to lie in the correct ball park. For a final conclusion concerning these 



matters, a reanalysis of presently available dimuon neutrino data is obviously mandatory! 



It is also interesting to check the above statements by an independent quantitative 
test sensitive to s(x, Q 2 ) such as for example the combination ^F£ N — SF^ which is given 
approximately by 

b -F» N {x,Q 2 ) - 3F? N (x,Q 2 ) ~ xs(x,Q 2 ) (6) 

where the charm contributions, the m 2 /Q 2 corrections and the NLO q- and g— induced 
contributions are rather small and, furthermore, \V CS \ 2 ~ 1 and \V c d\ 2 — 0. In fig. 3a we 



compare various LO and NLO results for xs(x, Q ) in eq. (6) with the published [12] and 



more recent but preliminary [jnj ^vN data and the NMC (deuteron) fiN data JT31. It 
should be kept in mind that the neutrino data refer to a Fe-target and are therefore very 
sensitive to nuclear (EMC) corrections in the small-x region: Only the preliminary (un- 
published) neutrino data ||13|| , which are larger than the published ones |12| in the small- x 



region, disagree with the approximate predictions, eq. (6), in fig. 3a. That this latter 
approximation is indeed sufficiently accurate is demonstrated in fig. 3b where xs(x, Q 2 ) 
is compared with the full NLO result for ^F% N — ZF^ which has to be calculated in the 
following way (note that F% always refers to an average over v and v). In the FFS we 
have 

F? N (x,Q 2 ) = x £ \(q' + q>)(x,Q 2 ) + [((q' + q')*C!)(x,Q 2 ) 
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with q' = |(1 + \V ud \ 2 )q+ \\V us \ 2 s, using \V ud \ 2 + \V US \ 2 = 1 with \V ud \ 2 = 0.9743, and 
where T% is given in eq. (2) with the replacement s' — ► |(s / + s'). The massless coefficient 
functions C|' 9 are standard, see e. g. ref. ||], and the convolutions are defined by 

(q*C)(x,Q 2 ) = J^q(z,Q 2 )c{^) . 

The well known expression for F% N [|J is appropriately modified for an isoscalar tar- 
get, with the charm contribution F 2 CC calculated according to the 7* g — * cc fusion process 



etc. [0] as described, for example, in |1. In the 'variable flavor' scheme |5|, where intrin- 
sic charm densities c(x, Q 2 ) are generated radiatively by the ordinary massless evolution 
equations, we have 

I ( 5 _ 3 f ,n\ = (fl _ c)(Xj g2) + _ c) + c|] (X) g2) _ (g) 



In view of the preliminary and contradicting nature of the nuclear-shadowing cor- 

(-) 

rected CCFR data for F 2 U used in fig. 3, a decision concerning the (dis) agreement with 
theoretical QCD predictions must obviously be postponed. According to our results in 
figs. 1 and 2, implying strongly that snlo is similar in size to s^o, and the ones in fig. 3 
which imply that the inclusion of the finite part of W + g — > cs and the corresponding 
photon induced 7*5 — > cc in conjunction with present NLO strange quark densities [4-6] 
do not change significantly the simple LO results, the theoretical predictions are rather 
constrained and unique. Furthermore, the results in fig. 3 again support our previous 



conclusions [0] concerning the perturbative stability H [Hd of the charm production rate 
as calculated in perturbative fixed order a s , i. e. in the FFS. A similar analysis was car- 
ried out in [16|| where different conclusions concerning the magnitude of the gluon induced 



contributions are presented: These results are almost a factor of two larger than the full 
NLO results at x = 10" 2 in fig. 3b since a factor of two error seems (due to the lack of 
explicit formulae in |16| it is not possible to trace its exact origin) to be present in the 
calculation of the W + g —> cs contribution. Therefore, if the enhanced preliminary uN 
data at x ~ 0.1 as shown in fig. 3 are confirmed, the discrepancy between these data 
and the rather solid and unique theoretical results, taking into account the rather well 
understood dimuon data as well, will constitute a major problem which cannot be solved 
within our present understanding of the so far successful perturbative QCD. 
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Appendix 



The fermionic NLO coefficient functions Hf for heavy quark (charm) production in eq. (2), 
calculated from the subprocess W + s — > gc, are given by 

12 i ,^,2 
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The coefficients in (A2) for i = 1,2,3 are given in Table 1 where a misprint in H con- 
cerning A 2 was corrected. 



Table 1. Coefficients for the expansion of hi in (A2) 
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The gluonic NLO coefficient functions Hf for heavy quark (charm) production in eq. (2), 
as calculated from the subprocess W + g — > as, are given by 



where P$(z) 



l ~ 3 

1 
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and 



hf(z, A) = C + C M z{\ -z) + C 2)i + (1 - A) z L x (C 3)i + A z C 4>i ) 



(A6) 



with 



C = P${z) [2 ln(l - z) - ln(l - Az) - lnz] . (A7) 
The coefficients Ck,i are given in Table 2 and differ from those in || where the older 



convention |L7j has been adopted of counting the gluonic helicity states in D = 4 rather 
than in D = 4 + 2 e dimensions. The latter convention |TE|] is the one chosen to define all 
modern NLO parton distributions. 

Table 2. Coefficients for the expansion of hf in (A6) 
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6A 
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Note that in the limit A — > 1 (m c — ► 0) the iff' 9 reduce, apart from the obvious 
collinear logs, to the massless MS coefficient functions Cf g H, [L8| . 
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Figure Captions 



Fig. 1 LO and NLO predictions for £s e // defined in eq. (4), using the GRV [|J and 
CTEQ3 || parton densities. The dotted curves refer to using just the NLO strange 
quark contribution in eq. (2) with all 0(a s ) terms neglected. Thus the differences 
between the dashed and dotted curves illustrate the differences between the LO 
and NLO strange sea densities, respectively. The values of Q 2 vary between 2.4 
to 43.9 GeV 2 according to the experimental averages || for 0.015 < x < 0.35 and 
E v = 192 GeV § has been used. 

Fig. 2 NLO results using the NLO strange sea density of CCFR 0. The subtraction 
term (SUB) is defined in (5) and an acceptance correction factor A = 0.6 has 
been used [llj. The analysis was performed with the original subroutines/matrix 
elements of CCFR ||11|| ; if the charged current structure functions of GGR f|| are 
used instead, the results are similar. The shaded area refers to the CCFR 'data' 
[§], calculated according to (4'), where the CKM suppressed contribution in (1) has 
been subtracted from the measured full cross section by assuming specific up and 
down quark densities 0, [□]]. The dashed curve corresponds to the original CCFR 



fit analysis 0, [Li 

Fig. 3 (a): LO and NLO GRV § and MRS (A) § predictions for xs(x,Q 2 ) which ap- 
proximates \F% N - 3F£ N in eq. (6). (b): Full NLO result for $F$ N - 3F£ N in 
the FFS (dashed dotted curve) using eq. (7), and the short-dashed curve shows the 
corresponding result with the W + g — > cs contribution turned off. The factorization 
scale chosen is fi 2 = Q 2 + m 2 . The solid curve for xs, being the same as in (a), is 
shown for comparison. The full NLO MRS (A) result in the 'variable flavor' scheme 
is based on eq. (8). Both CCFR |T2| (circles) and preliminary CCFR |T3| (squares) 
uN (Fe-target) data are corrected for nuclear shadowing effects, whereas the NMC 
fiN data |14| have been obtained from a deuterium target. 
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Fig. 3a 
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Fig. 3b 



